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A THEORETICAL STUDY OF THE ACOUSTIC IMPEDANCE OF ORIFICES 


IN THE PRESENCE OF A STEADY GRAZING FLOW 
by Edward J. Rice 

ABSTRACT 

An analysis of the oscillatory fluid flow in the vicinity of a circular orifice 
with a steady grazing flow is presented. The study is similar to that of Hersh 
and Rogers but with the addition of the grazing flow. Starting from the momen- 
tum and continuity equations, a considerably simplified system of partial dif- 
ferential equations is developed with the assumption that the flow can be de- 
scribed by an oscillatory motion superimposed upon the known steady flow. 

The equations are seen to be linear in the region where the grazing flow ef- 
fects are dominant, and a solution and the resulting orifice impedance are pre- 
sented for this region. The nonlinearity appears to be unimportant for the 
usual conditions found in aircraft noise suppressors. Some preliminary con- 
clusions of the study are that orifice resistance is directly proportional to 
grazing flow velocity (known previously from experimental data) and that the 
orifice inductiv^e (mass reactance) end correction is not a function of grazing 
flow. This latter conclusion is contrary to the widely held notion that grazing 
flow removes the effect of the orifice inductive end correction. This conclu- 
sion also implies that the experimentally observed total inductance reduction 
with grazing flow might be in the flow within the orifice rather than in the end 
correction. 

INTRODUCTION 

The acoustic impedance of perforated materials in the presence of a 
steady grazing flow is an essential component of noise suppressor research. 
This impedance is used as a boundary condition in all duct sound propagation 
studies. 

Several investigations of the grazing flow impedance phenomenon have 
been performed (Refs. 1 to 11) that lead to impedance models being generated 
for use in impedance calculations. Theoretical studies of a basic nature have 
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been performed (Refs. 12 and 13) but mainly without grazing flow. Recent 
work has been performed which is exposing the basic physics involved in the 
resistive component of the grazing flow impedance (Refs . 8 and 14) , however, 
the physics of the imaginary or reactive component are yet unclear. 

A valuable study has been performed (Ref. 15) which considered the 
acoustic impedance, both linear and nonlinear, of an orifice without grazing 
flow. This study is of particular interest since it was derived directly from 
the Navier-Stokes equations and was thus very basic in nature. Closed form 
solutions (exact for the linear and approximate for the nonlinear) were ob- 
tained for both the real and imaginary component of the acoustic impedance. 

The work reported in this paper represents an extension of Reference 15 
in that the steady grazing flow is added. A spherically symmetrica] flow per- 
turbation is assumed to be impressed upon the known unsymmetrical steady 
flow. Where the grazing flow effects dominate, the equations are found to be 
linear thus }delding closed form solutions (with suitable simplifying assump- 
tions) for acoustic pressure and velocity. The expression for acoustic impe- 
dance is also presented and the resistive part is compared to acoustic resis- 
tance models derived from experiments. The steady orifice flow solutions 
are also obtained and calculated pressures in the vicinity of the orifice are 
compared with measured pressures. 


SYMBOLS 

2 

Anl area dominated by nonlinearity, m 

O 

Atot area, m“ 

c speed of sound, m/sec 

d orifice diameter, m 

F(t) function of time (see Eq. (30)) 

f(t) function of time (see Eq (32)) 



grazing flow Mach number 
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Mq orifice Mach number 

^ exciting pressure amplitude far from orifice, N/m^ 

p nondimens ional perturbation pressure, p'/^ 

p* complete dimensional pressure (see Eq. (5)), N/m 

O 

p average or ambient pressure, N/m 

2 

p’ perturbation in pressure, N/m 

p nondimensional static perturbation pressure 

s 

o 

p* dimension form of p , N/m 

s s 

p* same as ^ except for steady orifice flow, N/m^ 

dynamic pressure % 'p V^), N/m^ 
r nondimensional radius (r*/d) 

r* dimensional radius, m 

Tq r at orifice radius (rQ - ^ ) 

t nondimensional time (ait*) 

t* time, sec 

Uq magnitude of radial perturbation velocity near orifice, m/ sec 

u nondimensional radial perturbation velocity (u'/Uq) 

u* complete dimensional radial velocity (see Eq, (6)), m/sec 

u radial component of grazing flow velocity (see Eq. (1)), m/sec 

u' radial component of perturbation velocity, m/sec 

Ug same as u but for steady onfice flow 

u^rf orifice perturbation velocity based on flow rate and orifice area, m/sec 

u|^^ orifice perturbation velocity at vena-contracta, m/sec 

grazing flow velocity, m/sec 
V 9 c’omponent of grazing flow velocity, m/ sec 
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V* complete dimensional 0-component of velocity (see Eq. (7)), m/sec 
v' 0-component of perturbation velocity, m/sec 

w (^-component of grazing flow velocity, m/sec 

w* complete dimensional </>-component of velocity (see Eq. (8)), m/sec 
w* <jP-component of perturbation velocity, m/sec 

5 phase angle between pressure and velocity 

^ dimensionless acoustic impedance 

0 polar angular coordinate (see Fig. 1) 

X wavelength, m 

Q 

p* density, kg/m 

p average uniform density, kg/m 

q 

p' density perturbation, kg/m 

cp azimuthal angular coordinate (see Fig. 1) 

oj circular frequency, rad/sec 


THEORETICAL DEVELOPMENT 


In this section first the differential equations will be presented in dimen- 
sionless form. Next, they will be simplified by a magnitude analysis and finally 
the solutions will be derived. 


Differential Equations 


The geometry of the system considered here is shown in Figure 1. Both 
rectangular and spherical coordinate systems are shown centered at the ori- 
fice center. The steady flow velocity is shown parallel to the x axis and al- 
though uniform in the rectangular coordinates, in spherical coordinates the ex- 
pressions are; 
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u =» sin 6 cos <p 

(1) 

u = -V_ cos d cos (p 

(2) 

w = sin 

(3) 


The differential equations, for the purpose of brevity, will be given with 
all of the assumptions inserted. The variables will be considered to be given 
by the sum of a steady component and a time var 3 'ing component as, 


p* = p 

+ p' 

(4) 

P* = P 

+ p’ 

(5) 

u* = u 

+ u' 

(6) 

l> 

II 

» 

> 

+ V' 

(!) 

w* = w 

+ w' 

(8) 


Asterisks are used here to denote complete dimensional quantities for the 
d^endent variables and dimensional quantities for the coordinates. The 
isentropic relationship is assumed valid so that 

p' = c^p' (9) 

The radial coordinate and time will be nondimensionalized as follows, 


r* - rd 

(10) 

CO 

(11) 


The time varying pressure is normalized by the peak far-field pressure 



6 


and thus using Equation (9) 


P' 



(13) 


The time varying radial velocity is normalized by the amplitude near the 
orifice (Uq) thus 


u' = UqU (14) 

The perturbation velocity (u'), the complete velocity (u*), and two components 
of the steady grazing flow (u, v) are shown in the two-dimensional sketch of 
Figure 2 . The simplest possible solution will be sought for the velocity and then 
tested to see if it has provided any information or insight into the physics of the 
grazing flow impedance. With this in mind, only axisymmetric perturbations 
of the velocity will be considered and thus, 

v' = w' - 0 (15) 

It should be noted later in the development of the equations that even though the 
velocity perturbation is axisymmetric, the pressure distribution around the ori- 
fice is asymmetric The viscous terms in the momentum equations will not be 
considered. 

With the above assumptions, the equations of motion can be written as 
follows; 


Co» tinuity 


P\ I o)d\^ 





1 8 


2 8r 


r 


(r^u) 





V ^ ^ ^ 1 9p 

c r 80 c r sin 0 dcp 


- 0 


(16) 
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r-Momentum 



(y>-Momentum 


1 ^,^0 

pencil ( (P \ d(P c 

' ' (1 +-^ Pjsin e 


(17) 


(18) 


(19) 


Although considerably simplified from the complete equations, Equations 
(16) to (19) are still much too complex for a simple closed form solution. An 
order of magnitude analysis on the coefficients must now be made in a manner 
similar to Reference 15 except that must be considered. Only the coeffi- 
cients need to be considen d since the variables themselves are of order unity 
due to the nondimensionalization. From several references (Refs. 4, 6, 7, 8, 
10, and 14 e. g.) the grazing flow resistance can be approximated by 


PcUq c 


( 20 ) 


Thus 



= « 1 


( 21 ) 



since is at most 0.6 or 0.7 for the usual acoustic liner and Mq is much 
less than this. Also 


isid = j ^22 

c X 

provided the wavelength (\) is much larger than the orifice diameter (d) which 
is almost always the case in a practical acoustic liner. The first term in the 
continuity Equation (16) is probably the smallest of all with the third term be- 
ing next smallest in magnitude (essentially Mach numbers to the third power 
when the terms in the square brackets are considered) . Also since p is of 
order unity the term ^Pp/'pc^ can be dropped where it is added to unity. Thus 
the problem to be solved is shown to be incompressible and the equations have 
been reduced to, 


— (r^u) = 0 (23) 

9r 


yc J 9t I c c j dr \ocUqJ dr 
I PcUq J 96 c 
I PcUq j sin 6 9yJ c 


(24) 


(25) 


(26) 


Note that Equation (24) is nonlinear in u and we will restrict the solution to the 
region in which the grazing flow velocity dominates the orifice velocity. An 
approximate nonlinear solution can probably be found following the method of 
Reference 15 but that will not be attempted here. A discussion of the range 
of validity of the linear solution will bo expanded upon in a later section. With 
this restriction in mind, Equation (24) can be written as 
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The boundary conditions to be applied involve the excitation pressure im- 
posed upon the system in the far field and tiie inevitable separation of the flow 
at the upstream edge of the orifice for inflow. This separation implies that 
for sufficiently short orifices (no reattachment) the ambient or back cavity 
pressure will be felt at die upstream orifice edge. This separation region is 
shown in Figure 2. For outflow the orifice boundary condition would be modi- 
fied and it is thou^t that the solutions could be made to model outflow. This 
is not done in this paper and the solutions which follow are intended for inflow 
only at this time. The boundary conditions can be expressed as follows 

p - e^^ as r-*- 00 (28) 

p = 0 at r=rn, 6 = <p = 0 (29) 

" 2 


Solution to Differential Equations 


The solution to the differential equations follow in a similar manner to that 
of Reference 15. Equation (23; can be immediately integrated to give, 



(30) 


where the negative sign was chosen since the velocity is opposite to the radial 
coordinate for inflow. The function F(t) d^ends upon t only since u was 
assumed not to be a function of 6 or <p. Using Equation (30), Equation (27) 
can be written, 


9r 






(31) 


Equation (31) can be integrated to obtain, 

(iiA A ^ _ u A F + 

\c Jr dt c j,2 \^pcUyy 


( 32 ) 


P "" f(t) 
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Use Equation (28) in Equation (32) as r-«- « to obtain, 


f(t) = 



(33) 


where f(t) can be only a function of time due to the boundary condition of Equa- 
tion (28). Thus, Equation (32) can be written. 


A ^ 
^c y r dt 


H A F + / (p - e^s = 0 

c 1 


(34) 


In solving for F(t) we will seek solutions of the form, 

F(t) = (35) 

Substituting Equation (35) into Equation (34) and imposing the upstream boundary 
condition (Eq. (29)) at the orifice edge, after considerable arithmetic, A and 6 
can be determined such that 


and 



tan 5 


wdr 


0 


V 

OO 


(36) 


(37) 


Using Equations (35), (36), and (37) in Equation (30) yields, 


- 7(V - iwdr„) 

QO U 


pU, 


v2 t (o,dr(,)2 



e 


it 


u = 


r 


(38) 
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and tiien Equation (34) will give, 



(39) 


It should be noted here that the velocity and pressure Equations (38) and (39) 
also satisfy the 0 and <p momentum Equations (25) and (26). 


ACOUSTIC IMPEDANCE 


Before calculating acoustic impedance one more point must be made. 
Usually Impedance is based upon the velocity calculated from orifice flow rates 
and orifice area. Note that the area of the hemisphere at r = rQ is twice the 
orifice area. The value of u calculated at r = rQ will thus be doubled when 
used to calculate impedance. 

The acoustic impedance w.Ul be defined as the far field perturbation pres- 
sure divided by the orifice velocity and then normalized by pc or 




_p> - 

P^“orf (2pcUou)j,^ 


(40) 


Using Equation (38) evaluated at r = rQ = l/2 in Equation (40) yields, 




( 41 ) 
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The negative sign occurs in Equation (41) because the velocity for inflow is 
in the negative r direction with the coordinate system centered on the orifice. 

A dicsusslon of Equation (41) along with a comparison with esqperimental 
data will be made in a later section. 


Steady Orifice Inflow Theory 

The theory of the preceding sections can be easily converted to a steady 
orifice inflow theory. This is useful since the theory can then be compared 
to the more easily obtained and thus more numerous steady flow data. This 
conversion to steady flow can most easily be obtained by letting oj -*■ 0 (thus 
also eliminating e^^ = e^^*) in Equations (38) and (39). These equations 
thus become, 




and 


Pg = 1 - ( — j sin 6 cos ip 


(43) 


Some interesting result" can be derived from the equations. First notice in 
Equation (43) that at the upstream edge of the orifice (r ~ r^. 6 - ir/2, (/? = 0) thai 
p = 0 which of course was the boundary condition reflecting flow separation at 
this point. The back cavity pressure is felt here. However at the downstream 
edge (r - Tq, 0 ■■ 7r/2. “ tt). Pg = 2 or twice the far field pressure, In fact if 

the pressure is integrated over a hemisphere arourd the orifice to find an average 
pressure, it is fo'ond that this average is equal to th»" far field driving pressure. 
The pressure has just redistributed due to the grazing flow. The pressure in the 
high pressure region of the flow must yet be relieved in passing to the back cavity 
which requires a further expansion. In olher words, the average pressure at 
r = rQ still has a magnitude The fluid must expand to have zero pressure as 
it passes to the back cavity side of the orifice. The velocity will thus increase 
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to a maximum at ttie vena-contiacta. Hiis will be analyzed using a one- 
dimensional approach to find an average vena-contracta velocity, u' , and 
then a disdiarge coefticient will be calculated. 

The fuitiier expansion in the orifice can be estimated from the Bernoulli 
equation. 





(44) 


where UqU^ is the dim^isitmal perturbation velocity at r = r^, and the pres- 
sure diY^ is from to zero in the back cavity. Using Equation (42) d ere 
results. 


Fr^m flow continuity considerations (use Eq. (42)) to calculate flow rate through 
die hemispherical area) the full flowing orifice velocity (based on orifice area 
and flow rate) is. 


U' 

orf pV 


(46) 


The orifice discharge coefficient is defined here as t ie ratio of the full flowing 
orifice velocity to the vena-contracta velocity and c; n then be calculated from 



where ^/2pV^ 


was assumed to be small compared to unity. 
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This result has been found empirically in Reference 14 and will be dis- 
cussed further in the next section along with pressure measurement compari- 
sons using Equation (43) . 


DISCUSSION OF RESULTS 


In this section the previously derived theory will be compared to data 
where it is available. The acoustic resistance is well established by experi- 
mental data. The acoustic inductance is less defined since measurements 
lump both the orifice end effects and within the orifice ma&s effects into one 
inductive quantity without any way of separating them. Also the steady orifice 
inflow theory will be compared to experiments. 


Acoustic Impedance 

Acoustic resistance . - The real part of Equation (41) (resistance) com- 
pares well with published experimental results. It is well established that 
resistance depends upon the grazing flow Mach number. 

The coefficient of in Equation (41) is 0. 5 which compares favorably 
with 0. 3 from References 4 and 7 for arrays of orifices and with 0.7 from 
Reference 14 for a single orifice. It is not intended that these theoretical 
results replace the empiricisms since there is much more woric to be done 
on the .letails of the flow within the orifice, but only to show that the theory 
appears reasonable. 

Acoustic reactance . - The theoretical results for orifice reactance are 
not so easily compared to data as was resistance. Since the major conclu- 
sions of this paper concern the reactance, it is best to first put the theory 
and the actual flow into proper perspective. Figure 2 shows sketches of the 
inflow to the orifice for conditions such that the grazing flow is completely 
dominant as assumed in the present theory. The flow patterns of figure 2 
have been established for steady orifice flow in Reference 14 and for oscil- 
latory orifice flow in Reference i6. The hemisphere of radius equal to 
the orifice radius is shown on Figure 2 and the solutions given in this paper 
are expected to apply outside of this hemisphere. The total inductance is 
determined by adding the inductance provided by the fluid within this hemi- 
sphere and within the orifice itself to the inductive orifice end correction 
given by Equation (41). Going back for a moment to zero grazing flow and low 
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pressure amplitudes, the orifice inductive end correction is given by 0.85d 
(Ref. 18) which would give 0. 425d at each end of the orifice. Equation (41) 
accounts fc over half (0. 25<b of this correction which is outside of the 
hemisphere shown in Figure 2 with the remainder neing within this hemi- 
sphere. Hard conclusions can thus be made about only half of this mass, but 
it is suspected that the conclusions will apply to the bulk of the attached mass. 

Experimental data (Refs. 6, 8, and 17) show that the total inductance is re- 
duced with a grazing flow. Since Equation (41) shows the inductive orifice end 
correction to be independent of grazing flow velocity, then the reduction in in- 
ductance which is observed must be coming from the fluid within or very close 
to the orifice. This conclusion is contrasted to the idea in which an orifice 
slug flow exists with an attached mass end correction which is blown away by 
the grazing flow. 

Additional study must be devoted to the flow in the orifice. The orifice 
flow (with a grazing flow velocity) is Mt a slug which accommodates more mass 
flow by an acceleration of the slug. Instead the flow may have very nearly a 
constant velocity (see Ref. 16) with additional mass flow accommodated by an 
increased flow area. A sketch of these flow patterns is shown in Figure 3. It 
may be that a constant velocity, variable area type of flow has essentially no 
inductance at all. This latter point is currently being studied. 


Steady Orifice Fiv.w 

For steady inflow into an orifice some discharge coefficient data are 
given in Reference 14. At low orifice velocity to grazing flow velocity ratios 
a fit to that data can be made as 



This result compares favorably with the theoretical result of Equation (47). 

In Reference 14 some pressure tap data was given in the vicinity of the 
orifice for both upstream and downstream locations. The data is given as 
(p* - p*)/q «3 For upstream pro.ssure. Equation (43) gives 
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P 


s 


= 1 



(49) 


where 9 = it/2, <p = 0 were used. 

Equation (49) can be converted to (using Eq. (46)), 



(50) 


This is plotted in Figure 4 along with the results from Refer^ce 14. The 
radius values used were r/r^ = 1. 47 which was determined by fee dimensions 
shown on the figure insert. The agreement between the theory and fee data are 
seen to be reasonably good. 

For downstream pressure taps fee analog of Equation (50) is, 



ipv^ 

2 ~ 


where $ -- jr/2 and <p = i: was used in Equation (43). This result is also shown 
in Figure 4 along with data from Reference 14. The radius values here were 
r/rQ=1.20 determined from the dimensions from the insert on Figure 4. The 
odd behavior of the data at very low orifice flow was caused by the vortex stand- 
ing in the orifice, which raised the streamlines above fee orifice, and slightly 
depressed the downstream pressure. Once this depression effect is overcome 
by going to larger values of u^^^ the theory and the data agree quite well. In 
fact, the agreement for both curves in Figure 4 extends beyond fee expected 
range of validity of the solution which will be discussed next. 


RANGE OF VALIDITY OF THE LINEAR SOLUTION 


Recall that the r-momentum Equation (24) was actually nonlinear and that 
this nonlinearity was ignored in the solutions of this paper. The coefficient of 
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the second term in Equation (24) may be written, 

(u + UqU) = - sin 0 cos + UqU (52) 

In the area where the first term of the equation dominates the second term, the 
present solution might be expected to approximate the complete solution. This 
may be true provided the area of influence of the nonlinearity is small. The 
coordinates of the flow area at which the two influences are equal is given by. 


UqU 

sin d cos <p = 


(53) 


The nonlinear solution is approximated from the results of Reference 15 (and 
ignoring die time dependence) as 


UqU 



(54) 


thus Equation (53) becomes 


sin 6 cos <p fs 



(55) 


The ratio of the flow area dominated by the nonlinear solution to that of the en- 
tire area may be estimated by (assuming is very small), 

00 


^NL ^4 j~^ 

^TOT pV^ 


(56) 


Mass flow rates in the two regions (rather than areas) could also have been used 
to indicate the importance of the nonlinearity. Using Equation (46), Equation (56) 
becomes, 
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‘NL 


'TOT 



(57) 


Thus it is presumed that as long as the quantity on the ri^t side of Equa- 
tion (57) is small, then the nonlinearity can be ignored. Note from Equation (52) 
(when sin 6 cos « 0) that the nonlinear term should dominate when 0 = 0 or 
<p - ff/2 or 3ir/2 which is above and to the sides of the orifice (in the y-z plane 
of Fig. 1) . Also note from Equation (43) (again sin 9 cos » 0) that this is pre- 
cisely the region in which the linear solution would predict that the pressure is 
equal to the far-field driving pressure (both equal to one) . Thus with an apparent 
lack of driving pressure difference in the nonlinearly dominated area, the linear 
solution may be valid well beyond the criterion given above (the ri^t side of 
Eq. (57) being small) . This may explain the agreement between the linear theory 
and the data shown in Figure ' well beyond the range at first expected. 

Thus it is anticipated that for the usual acoustic pressures and grazing flow veloc- 
ities found in aircraft suppressors, the nonlinearity and thus the nonlinear acous- 
tic resistance due to the acoustic pressure will be ne^igible. 


CONCLUDING REMARKS 


The main conclusion of this paper is that the orifice inductive and correction 
or attached mass is not a function of grazing flow velocity. This was shown by 
calculating the acoustic impedance in Equation (41) and observing the inductance 
term. This conclusion could also be drawn from the governing differential Equa- 
tion (27). Only the first term in this differential equation can contribute to veloc- 
ity components which are out of phase with the pressure and this term does not 
contain the grazing flow velocity. This approach (using (he differential equation) 
is more general since now also the orifice outflow can be included if it is as- 
sumed the small separated flow leaving the orifice near the wall on the down- 
stream side does not appreciably affect the reactive part of the solution. 

When the nonlinearity (small for Uqj./^oo than unity) is restored 

to the differential Equation (24) then the conclusions must be modified somewhat. 
The nonlinearity can produce reactive changes as shown in Reference 15 but 
these effects, especially their dependence upon grazing flow, will be unknown 
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until the nonlinear solution is obtained. However, it has been argued here that 
this nonlinearity is probably not important except at extremely high pressure 
amplitudes and small grazing flows for which would become large. 

If the inductance change with grazing flow does not come from the attached 
mass (orifice inductive end effect) then it must be coming from the flow within 
the orifice itself. This region is a prime candidate for this effect since the 
flow within the orifice may have close to a uniform velocity with mass flow 
changes occurring due to a variable flow area. 

The acoustic resistance was calculated using the theoretical results pre- 
sented here. The resistance was shown to depend linearly upon the grazing 
flow velocity with a magnitude agreeing quite well with published data. Since 
the nonlinearity in the differential equation does not seem to play an important 
role, this implies that the resistance will be dominated by the grazing flow 
and will not be sensitive to the acoustic pressure. 

It should be recalled that the solutions to the differential equations pro- 
vided here are for the flow perturbations. To obtain the complete pressure 
and velocity solutions, these perturbations must be added to the average pres- 
sure and uniform grazing flow velocities. 

The complete solution to the orifice impedance problem is yet to be ob- 
tained. This will involve the coupling of the solutions on the cavity side of the 
orifice, the complicated three-dimensional orifice flow, and the grazing flow 
side of the orifice for both inflow and outQow. 
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Figure 1. - Geometry used In orifice flow analysis. 
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Figure 2. - Sketch of orifice inflow with a grazing flow. 
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Figure 3. - Sketch of orifice flow patterns. 



Figure 4. - Comparison of measured and theoretical static pressures near 
an orifice. 



